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00 ! Abstract 

We deal with a system of two coupled differential equations, describing the 
X ', evolution of a first order phase transition. In particular, we have two non- 

^ I linear parabolic equations: the first one is deduced from a balance law for 

entropy and it describes the evolution of the absolute temperature; the other 
one is an equilibrium equation for microforces and it regulates the behaviour 
of a scalar phase parameter. Moreover, the second equation shows a time- 
relaxation coefficient related to the time-derivative of the phase parameter. 

We prove well-posedness of solution to the given system, using a standard 

method of approximating problems; afterwards, we study the behaviour of the 

1/-^ ' system as the time-relaxation coefficient tends to zero: as a result, we find out 

op ■ that the original problem converges to a new problem, with a stationary phase 

equation. 



1 Introduction 



This paper deals with a system of two coupled differential equations, describing 
k> ! the evolution of a first order phase transition. Both equations are non-linear and 

H I parabolic: the first one is deduced from a balance law for entropy and it describes the 

evolution of the absolute temperature '&; the other one is an equilibrium equation for 
microforces, which are responsible for the phase transition process, and it regulates 
the behavior of a scalar phase parameter x- 
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The system is regarded in the parabohc cyhnder Q := Qx (0, T), with T > and 
n C R^, smooth and bounded, and it has following expression: 

dt{log^ + x) - M = g 

^JiXt-^X + i + ^\x)='^ and ied[5{x)- 

where g is a, thermal source, (9/3 is the subdifferential of an appropriate proper, convex 
and lower semi-continuous function /3 ((9/3 is indeed a maximal monotone graph) and 
o" is a C^ function with Lipschitz continuous derivative. 

The sum 8/3 + a' introduces a non-linear term related to some physical constraints 
(thermodynamically consistent) which can eventually occur in applications. 

Finally, the coefficient n represents a time-relaxation parameter and it is usually 
a small quantity in applications, compared with other physical quantities which act 
during phase-transitions. 

The PDE system is completed with initial and boundary data: a non- homogeneous 
Robin condition is provided for the absolute temperature t9 and a standard Neumann 
condition is set for the order parameter x- 

d^d + a'd = h on(91]x(0,T) 

(9^X = on(91]x(0,T) (1.2) 

log'i9(0) = log-i^o and x(0) = Xo in ^5 

where a is an arbitrary positive function, with appropriate regularity, and h denotes 
the boundary thermal supply; dy indicates the outward normal derivative. The initial 
data are given by two function df) and Xo defined on Vt and in particular we set -i^o > 
on Vl. 

The main difficulties in dealing with this system are the non-linear terms, in 
particular the logarithmic term and the maximal monotone graph (9/3. However, the 
presence of the logarithm enable us to conclude straightforward the positiveness of 
the temperature d, once the existence of solution has been proved in some suitable 
functional space. 

To obtain the result of existence of solutions, we have to go through a double 
approximation of the problem: first we regularize the non-linear terms and the ini- 
tial data and subsequently we proceed with a Faedo-Galerkin procedure to solve 
the approximated system. Finally we attain to the existence by using appropriate 
compactness results. 

After that, we are able to prove the uniqueness of solution to our problem. 

The second target of this paper is the study of the behavior of the system as the 
coefficient /x tends to 0. In particular, we will find out that the original problem (P^^), 



with /i > 0, converges to a problem (Pq), where the time derivative of x does not 
appear in the phase equation. 

The study of this asymptotic behavior has remarkable interest in many physical 
situations, where usually the fi parameter is smaller with respect to other physical 
quantities (such as the interfacial energy coefficient, which is related to the laplacian 
term in the phase equation). 

Since we have a stationary equation for the phase parameter, as /i \ 0, we 
loose time regularity properties for x and consequently for the logarithmic term, too. 
Thus, we need to impose further hypotheses on the non-linear terms P and a and 
to introduce a generalization of the function log'i? in the space {H^{Q))', in order to 
prove the existence of solution to the limit problem [Pq). 

This system was first introduced and studied by Bonetti et al. in [3] and |1]. 
In those articles the mathematical model of the system above (with extra terms 
involving a thermal memory kernel) is explicitly derived and existence and uniqueness 
of solution is proved with /i > fixed and with Dirichlet boundary condition for the 
variable d. Moreover, the latent heat is thought as an arbitrary function of X)iiot 
constant, as we assume in this paper. 

On the other hand, the use of Robin boundary condition in this paper introduce 
new difficulties in our system in proving the existence of solutions, since we have to 
deal with boundary terms which have to be studied using appropriate trace theorems. 

A similar model, which shows an entropy balance equation for the evolution of 
the temperature 'd, instead of the usual energy balance equation, has been studied by 
Bonetti, Colli and Fremond in [5J and by Bonetti and Fremond in [6]. In particular, 
both articles show explicit forms of non-linearities, while in the present article the 
non-linear term is represented by generic functions /3 and a, and the logarithmic term 
appears also under the laplacian; moreover, the article [5] shows a thermal memory 
kernel. 

Regarding the study of asymptotic behavior of the solutions as a certain coefficient 
tends to zero, we can quote the article of Gilardi and Rocca P], where, starting from 
the same model proposed in [1] and [3], the convergence to of the energy interface 
coefficient in front of Ax in the phase parameter equation is analyzed. 

We quote also the article of Colli et al. [9], where the asymptotic behavior as the 
time relaxation parameter /x \ is analyzed, but we point out that in this case a 
Caginalp phase-field model with memory is taken into exam; moreover, the equation 
regarding the behavior of the temperature d is derived from an energy balance (and 
not from an entropy balance as in this article) and does not show a logarithmic term. 

The paper is organized as follows: in Section [2] we provide a description of the 
mathematical model which is later studied in Section HI where the well-posedness of 



the problem is proved. Finally, in the last Section |5]we study the asymptotic behavior 
of the problem as the time relaxation parameter /i tends to 0. All the results of the 
present paper are stated in Section [3l 

2 The model 

Consider a first-order phase transition occurring in a smooth domain of the physical 
space fi C R^. The unknown variables of our problem are the absolute temperature 
f? € (0, +oo) and the order parameter x- Iii particular, x represents a local concen- 
tration or the rescaled proportion of one phase with respect to the other and it is 
related to the microscopical movements of particles. 

In order to grant thermomechanical consistency to the model, we need to intro- 
duce some constraints on x, which will be derived from suitable functions fi and a, 
as we will see later. 

The derivation of the model is mainly taken from the article by Bonetti et al. 
([3]), which is based on an approach independently proposed by Gurtin ([?]) and 
Fremond (|22])- We repeat it here for the reader's convenience. 

Assuming that macroscopic deformations do not occur during the process and 
that microscopic accelerations are negligible, we start stating the following balance 
law for microforces (compare with [22j): 

divH + h = B mVL . . 

H-z/ = o on an, ^ ' 

where H and B are interior microscopic forces and h is an external force acting on the 
body at microscopic level; in particular, we are assuming that no external contact 
force acts on the boundary of our domain Vt: this is the reason why we imposed a no- 
flux boundary condition for H. We can think of all these microforces as mechanically 
induced heat sources, which have to be taken into account while dealing with the 
first law of thermodynamics (in a local form): 

et = -divq + r + div(Hxt) + bxt = -divq + r + H ■ Vxt + Bxt, (2.2) 

with e the internal energy of the system, q the heat flux and r an additional heat 
source. 

At this point, we aim to study the thermomechanical consistency of the model and 
to find explicit expressions for the microscopic forces B and H. First of all, let recall 
the second law of thermodynamics in the form of the Clausius-Duhem inequality 

Vt > -divQ + g, (2.3) 



with 7] the entropy of the system, Q := c[/i) the entropy flux and g := r/i) the 
external entropy supply; we introduce also the free-energy functional \E' (see [5], [2T] . 
dUeH) 

vl/(t9,x) = ^(t9,x) + ||Vx|% (2.4) 

with ^ the density of the free energy for pure phases; |VxP takes into account 
local interactions between phases, while the constant £ > is the energy interface 
coefficient (see PHi Sec. 3 and 4] and [27, Sec. 3] for a thorough discussion on the 
argument). 

Using the Helmholtz's relation 

" = -IM • (2.5) 

and assuming that the entropy has the following expression 

r] = c,{l + \ogi^)+ix, (2.6) 

where Cg G R+ is the specific heat and i G R+ is the latent heat associated with the 
phase transition, we have 

*(^?, X) = -c.^log^ - Mx +[/3 + a] (x) + | |Vx|' , (2.7) 

with /3 : R — )■ [0, +oo] a proper, convex and lower semi-continuous function and 
(T : R — )■ R another (suitably regular) function. The functions we have right now 
introduced represent arbitrary constraints on the order parameter x; by suitably 
choosing /3 and a, the model may describe different types of phase transitions. 

Examples of possible non-linearities. The non-linear term P + a may have 
different significant expressions. 

If we are assuming that the two phases may coexist with different proportions (i.e. 
we are assuming the possibility of the so called "mushy regions"), we can reasonably 
impose the following constraint 

XG[0,1] (2.8) 

setting I — X equal to the proportion of the other phase; in particular, the values 
X = 0, 1 correspond to the pure phases, while for the in-between values both phases 
coexist at each point of the body. 

To force the order parameter to take only the required values, we can set /3 = /[o,i] 
the indicator function of the interval [0, 1], defined by /[o,i](x) = if x G [0, 1] and 



-^[0,1] (x) = +c>o elsewhere. The sub differential of (3, which will appear in the phase 
equation, is then the following maximal monotone graph 

^ < if X = 

^edIio,i]ix) ^^ <( e = ifxG(0,l) (2.9) 

e>o ifx = i- 

In this case, we can choose the following expression for the function a, which is quite 
common for solid-liquid phase transitions, 

a{x) = -iA + ^Mx{l-x), (2.10) 

where '&c is the critical temperature of the transition and M is the maximum value 
of a, attained at x = 1/2 (see pi Sec. 2.4]). 

Another important model is the so called "double-well" model, obtained by choos- 
ing the following form for the potential 

W = P + a = ^{x'-l)\ (2.11) 

which has two minima at the points x = =tl (pure phases) and a maximum at the 
point X = (transition point). 

In this case, the non linear term in the phase equation will be (3'{x) + ^'(x) = 
X^ -X- 

Combining the first and second law of thermodynamics ( (I2.2p and fl2.3j) ). we get 
the inequality 

^t = Ci - ^tT] - % 

< [-divq + r + H ■ Vxt + Bxt] - M - ^ -div (^|) + ^ 

= -^^r^ - iq • V^ + 5xt + H ■ Vxt, (2.12) 

which has to be identically fulfilled by any admissible process ^ = {^t,Xt, Vxt, Vi?). 
The first equality follows from the well known relation: e = \E' + i)r]. 

Let consider the processes ^ = (0, Xt, 0, 0) and ^ = (0, 0, Vxi, 0) and apply 
them to the equation fl2.12p : 

-T^Xt -Bxt<0 
f d^ \ (2.13) 



for each choice of Xt and Vxt- 

Under the hypothesis of small perturbations, we can assume the functionals B 
and H to be linearly depended on the dissipative variables: 

B = Bnd + BdXt and H = H„rf + i/^Vxi- (2.14) 

Combining these two expressions together with the previous inequalities (12.130 . 
we obtain the following state laws 

Bnd=^ = -M + ^ + a'ix), Bd>0, where ee9/3(x), 

d^ ^2.15) 

H«^ = a(V^ = ^^^' ^^^0- 

Moreover, we set 

Bd = fi>0 constant, H^ = 0. (2.16) 

Furthermore, recalling the general formula for heat conductors q = —k{'t})V'&, we 
set the heat conductivity A;(^) to be a linear function of the temperature: 

q = -fco^V^, with ko > 0; (2.17) 

this choice is quite standard for many dielectrics, like ice or water. 

At this point, we can rewrite the first law of thermodynamics (12.21) in the following 

way 

et = -divq + r + {Bnd + BdXt)Xt + H„d ■ Vxt 

= -Miv (^1) - Iq . V^ + r + BndXt + Bd |xt|' + ^nd ■ ^Xt- (2.18) 

On the other hand, thank's to relations f l2.15p . 



et = '^t + i^tV + ^r], 






+ dtr] + dr], 



= S„dXt + H„d ■ Vxt + %. (2.19) 

Then, by comparison, we have: "drit = — ^div (— fcoV'i?) — ^q ■ V'd + r + ft \xt\ ■ 

Now we are allowed to neglect the high order non linearities, thank's to the 
assumption of small perturbations, divide the equation by ^ > and explicit the 
entropy 77: 

dtics\og^ + ix)-koA^ = g, (2.20) 
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where g := r/-&. 

Finally, assuming that no external microscopic force acts on the body (6 = 
in equation (12. ip ). thank's to relations (I2.15P and (I2.16p . the microforces balance 
equations turns out to be the following parabolic non linear equation 

eAx = ~M + ^ + cx'{x)+f^Xt, (2.21) 

with e e dPix). 

The above equations are combined with suitable boundary and initial conditions. 
In particular, concerning boundary conditions and recalling (12. ip . we fix a homoge- 
neous Neumann condition for x'- 

li-u = ed^X = ondnx{0,T), (2.22) 

with z/ the outer normal derivative on the boundary of the domain; while, for the 
absolute temperature, we fix a Robin condition 

d^^ + a^ = h on 9fi X (0, T), (2.23) 

with a a suitable real-valued function defined on dQ and h a boundary heat source. 
The initial data are the following 

log^(0)=log^o and x(0) = Xo- (2.24) 

3 Main results 

In this section we show the rigorous formulation of our problem (ll.ip - (ll.2p and we 
state our main results. 

3.1 Well-posedness of problem (P^), with fi > 

Let consider a bounded, connected domain ^2 C R , with suitably regular boundary 
T := dn . Let T > be the final time and Q := n x (0,T) and S := L x (0,T) be 
the parabolic cylinder and the parabolic boundary, respectively. 
We set the following Hilbert triplet (see [T9] ) 

{V, H, V) = iH\n), L^n), iH\n)y) (3.1) 

with the usual standard norms and inner products. 



Let define the following linear, continuous operators 

B\V-^V\ {Bu,v)= Vu-Vvdx+ auvds yu,veV; 

■^^ ^^ (3.2) 

A:V-^V', {Au,v)= Vu-Vvdx \/u,veV; 

Jn 

with a G L°^(r): 

< a < a{x) < a a.e. on F, for some a, a G R+. (3.3) 

Remark 3.1. The operator B is symmetric and coercive, thus it defines a scalar 
product in V , which is equivalent to the standard one (thanks to Poincare's inequal- 
ity). From now on, we refer to || ■ ||y for the norm induced by B. 

The hypotheses on the source terms are the following ones 

geC\[0,T];H), 

heL\0,T;L^{T))r]W^^\0,T;L^{T)), h>0 a.e. in S; ^'' 

moreover, we define the operator w G C"^([0,T]; V) as 

{w{t),v) = / g{t)vdx+ I h{t)vds \/t G (0,T), \/ v G V. (3.5) 

in Jt 

The constraints on the phase parameter are given by the following functions 

/3 : R — > [0, +oo] is a proper, convex, lower semi-continuous function 

such that /3(0) = 0; (3.6) 

a G C^(R), a' is Lipschitz continuous, with Lipschitz constant cl- 

Hence, we consider the sub differential 9/3 of the function /3 which turns out to be a 
maximal monotone graph, with 9/3(0) 9 (see p^ Chap. II]). 

Remark 3.2. Thanks to Ii3. 6]) . we deduce the following polynomial growth for a 

k(r)| <c^(l + r2) \/reR. (3.7) 

Finally, the initial data has the following regularity properties 

^0 e L^{n), ?9o > a.e. in Vt and l/?9o e L°°(fi). ^ ' ^ 
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Moreover, these conditions imply < ^* < ^0(2;) < "&*, for a.e. x & fl, for some 

^„r eR+. 

In conclusion, our problem [Pfj,), with fi > fixed, turns out to be the following: 
we search for a pair (-(9^, x^) satisfying 

dt (log -df, + Xf.) + B^^ = w in \/' and a.e. in (0, T) (3.9) 

MXm + ^Xm + ^m + '^'(Xm) =^M a.e. in Q (3.10) 

e^ e 9/3(x^) a.e. in Q (3.11) 

log??^(0) = logt^o and Xm(0) = Xo,/. a.e. in fi. (3.12) 

Remark 3.3. In this formulation of the problem (P^), one can notice that some 
coefficients are missing, in particular the energy interface coefficient e, the latent 
heat i, the specific heat Cg and the heat conductivity k^. For the sake of simplicity, 
these physical constants had been normalized to 1. 

Moreover, thanks to the properties of the functionals A and 5, the following 
functional spaces are well-defined 



D{A] H) := A'^{H) = [u e H^{VL) | 9^m = } 

D{B;H) ■=B-\H) = {ueV\AueH, d^u + au = h} . 



Theorem 3.4. Let fi > and T > fixed. Under the assumption l[3.cl\) . (TO), ( I A*. 6]) . 



I{3.8\) . there exists a unique pair (-i?^, Xfi) C'lT'd there exists a selection ^^ such that 

^f,eL'^iO,T;V), ^^>0 a.e. inQ 

log^^ G L^{0,T;H)nH\0,T;V') 

X, e L'{0, T; D{A; H)) n H\0, T; H) ^ ' ' 

^, e L\Q) 

satisfying the problem h3.9\) - [37T^) . 

If the norms of all data, related to ^3.4^ and 1(3.3^) . are bounded by a positive 
constant M, then the solution {'&fj,,Xfj,j^fj.) satisfies the following estimate 

\\''^fJ'\\L^{0,T;V) "*" \\^^^''^fJ-\\ L°°(0,T;H)nm(0,T;V') '^ W^I^W L'2(0,T;D{A;H))nm{0,T;H) 

+ U,h.^O,T;H)<M' (3.15) 

where M' = M'{n,T,M). 

Moreover, the components {"t}^, x^) of solution continuously depend on data, in the 
following sense: if {gi,hi,'doi,Xo,fj.,i), "^ = 1,2, are two sets of data, whose norms are 

10 



bounded by two constants Mi and M2 respectively, then the corresponding solutions 
{^n^i, XtJ.,i) f'^^fi^ the following estimate 

jQt JQt 

+ / |V(Xm,i - Xm,2)|' + /i / IXm.iW - XM.2WI' + 111 * (^M,i - V)Wlly 
JQ Jn 

V0,fj.,l - V0,fi,2\\H + ll^'l ~ 92\\l2(^o,T-H) + 11^1 ~ ^2|lL2(o_T;L2(r)) (3.16) 



<M" 
yt G (0,T); «;/ierer7o,M,i := \og^oi+Xo,,^,i>fort = 1,2, antiM" = M"(fi,T, Mi, M2) G 

3.2 Asymptotic behavior as ^ \ 

Starting from problem (Pfj,), we can formally write down the limit problem (Pq), 
consisting of the first equation (13.91) integrated in time and of the second equation 
fl3.10p without the time derivative of the phase parameter x- 

Since we loose regularity estimates for the term dtX as /x \ 0, the logarithmic 
term logf? is less regular and takes values only in V'; thus, we have to modify the 
formulation of the problem and replace the term log ^ by a new unknown function ( G 
L'^{0,T]V') which generalizes the equality ( = logi?, through a suitable functional 
relation between ( and i). Basically, we define an ad hoc generalized logarithm as it 
has been done in 0, to which we refer for a detailed discussion on this argument. 

Definition 3.5. For ^ G L^{0,T;V) we define Log^ as the set of ( e L^{0,T;V') 
such that 

{C,e-^)+ [ ^Ij{^)< [ ^{6) \/eeL\o,T;r) (3.17) 

Jq Jq 

where ^(r) = r(logr — 1) if r > 0, ip{0) = and ^(r) = +oo if t < 0; we set 
D{Log) = {d e L'^iO,T;V) \Log^^iD}. 

At this point, we define "^ : L^(0,T; V) — )■ (— oo, +oo] as 



^i{v)= I ^{v) \/vEL'^iO,T;V). (3.18) 

Jq 



\E' is proper, convex and v G -D(^E') if and only if v is non-negative. 

Moreover, it is possible to prove that d"^ : L'^{0,T]V) — )■ L'^{0,T]V') is well- 
defined and Logv = d'^^v) (see I2j Section 4])- 



11 



The problem (Pq) can be written as follows 

Cit) + xit) + 1 * B'd{t) = I * w{t) + rio in l^', a.e. in (0,T) (3.19) 

C G Log?9 in V, a.e. in (0,T) (3.20) 

Ax + ^ + a'{x)=^ a.e. in Q (3.21) 

e e d(3{x) a.e. in Q, (3.22) 

being tjo := logf^o + Xo! the symbol * denotes the usual time convolution product 

{a*b){t):= a{s)b{t-s)ds, VtG[0,T]. (3.23) 

Jo 

To perform the asymptotic analysis as yU \ in problem (P^) (I3.9p - (l3.12p . we 
ask further regularity for the functional w (compare with (13. 5p ) 



geH\0,T-H) 

h e H\0, T; ^^(r)), h>0 a.e. in E, 

which imply w G H^{0,T; V). 

The hypotheses on initial data are the following ones 

XoeV, f3{xo)eL\Q) 

i3o e L~(fi), do>0 a.e. in fi e l/?9o e ^°°(fi). 



(3.24) 



(3.25) 



We assume also that the sequence of initial data {xo,At} related to problems (P^) 
satisfies 

Xo.m'^Xo in ^ and \\xoJ\v + \\l^iXo,f^)\\mn) ^ ^ V/x G (0, 1). (3.26) 

Regarding the functions (3 and a, we require the same properties as before (see 
(I3.6P ). together with some further hypotheses, which turn to be useful to prove 
the convergence of problem (P^) to (Pq). In particular, we deal with two distinct 
hypotheses. 

Hypothesis 1. Let assume a to be a linear function and f3 to satisfy the following 
condition 

/3(r) > cir^ - C2 for each r G D{(3), with ci, C2 G R+. (3.27) 
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Hypothesis 2. Let assume the property 

(6 + o^'(Xi)) - (6 + a'(X2))] (Xi - X2) > P (Xi - X2f , (3.28) 

for each Xi ^ D{df3) and for each ^j G 9/3 (xi), "^ = 1,2, for some positive constant p. 
Moreover, we ask that f3 + a has at least square growth 

/3{r) + a{r) > Civ^ - C2 for each r G D{/3), with Ci, Ca G R+. (3.29) 

We can notice that the second hypothesis is fulfilled if, for example, dP is strictly 
monotone and the growth of a' is dominated by d(3. 

Theorem 3.6. Let fi G (0,1). Assuming the same hypotheses of Theorem, \3.4\ let 
suppose that ^3.24 ), Ii3.25\) . Ii3.26\) hold and let assume either Hypothesis 1 or 2 on 
the functions /3 and a. 

Then, the solution to problem (P^) l^3.9\) - [37T^) . given by Theorem \3.4\ converges 
to a solution to problem (Pq) Ii3.19\) - f3.22\) . as /i \ 0, with respect to the natural 
topologies and with the following regularities 

i^eL\0,T;V) 

CGL-(0,T;n 

XeL\0,T;DiA;H)) 

^eL\0,T;H). 

Theorem 3.7. Under the same assumptions of Theorem \3.6i the components {'&,x) 
of the solution continuously depend on data; i.e. if {gi,hi,{}oi,Xoi), "^ = 1,2, are two 
sets of data, whose norms are bounded by constants Mi and M2 respectively, then 
the corresponding solutions {'di, Xi) fulfil the following estimates. Under Hypothesis 
1, we have 



(3.30) 



(Ci - C2)(^i - ^2) + / (6 - 6)(xi - X2) + II V(xi - X2) 



1 * (^1 - ^2){t)\\y < M ||r/o,l - T1qA\ + 11^1 - ^2||i2(o,T;//) + 11^1 



|2 



+ 



h II 

(3.31 



Vt G (0,T); where r]oi = hg^oi + Xoi, ^ = 1,2, and M = M(n,T,Mi,M2) G R^ 
while, under Hypothesis 2, we have 



<M 



(Ci - C2)(^i - ^2) + IIXi - X2lL2(o,T;y) + 111 * (^1 - ^2)(t)||^ 



2 ii2 ii2 

ho,l - ^70,211^ + \\9l - 92\\l2(0,T;H) + 11^1 ~ ^2 ||i2(o,T;L2(r)) 



(3.32) 



yte (0,r), with M = M{n,T,Mi,M2) G R^ 
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The following Sections H] and [5] are devoted to the proof of the previous results. 



4 Proof of Theorem 13.4 



In order to prove the well-posedness of (-P^), /i > 0, we proceed by a double approx- 
imation of the problem: first we introduce a family of regularized problems {{Pe)} 
depending on a positive parameter e; next, we prove the existence of a solution to the 
problem (Pe), e > fixed, applying a standard Faedo-Galerkin method. Afterwards, 
we let e tend to zero and deal with the original problem (P^). 

We base our proof on the guidelines of article [H Section 3, 4 and 5]. The main 
differences with respect to [1] are due to our boundary conditions: we will carefully 
detail this point in the proof. 

4.1 Hypotheses and preliminary results 

We consider the Yosida regularizations /3^ and log^ of the maximal monotone graph 
df3 and log respectively (see Reference [lEl Chap. II]), and define /3s, Log^ : R — )■ R 
by 

Ps{r) := / P's{s)ds and Log^(r) := £:r + log^(r). (4.1) 

Jo 

We can notice that both /3^ and Log^ are monotone and Lipschitz continuous. 

We need one more function, namely 



Is{r) := / sLog'^{s)ds, r G R, 
Jo 



(4.2) 



which is an approximation of the identity on (0, +oo). 

At this point, we state two lemmas which will be useful later. 

Lemma 4.1 (see Lemma 4.2, [1]). The function Log^ is differentiable, with derivative 
Lo(/^ such that 

e < Log'^ir) <e + -, WreR, (4.3) 

provided that e is small enough. 
Lemma 4.2 (see Lemma 4.3, [Ij). For all r ^ R the following inequality holds 

/.(r) < |r2 + 2r, (4.4) 

provided that e is sufficiently small. 
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Finally, we introduce regularized data {^oe} such that 

^oeeV \fe>0, < ^, < ^oe < r a.e. in Q, Ve > 0, 
■i^oe > "^0 in H and a.e. Q. 

Hence, the approximating problem (Pg) is the following one 



(4.5) 



dt (Log^'i?^ + Xe) + B^e = w in V', a.e. in (0, T) (4.6) 

fidtXe + Axe + f^'eiXe) + ^'(Xe) = ^e a.e. in Q (4.7) 

t9e(0) = ^0. and x.(0) = Xo,/. a.e. in n. (4.8) 



Theorem 4.3. Let e > fixed. Under the same hypotheses of Theorem \3.4\ and with 
the further hypotheses i4-5^ , there exists a unique pair [df^^iXs) such that 

^eeL\0,T;V)nH\0,T;H) 
XeeL\0,T-D{A-H))nH\0,T;H) ^ ' ^ 



and fulfilling the equations ([^.6])-([7^ of problem (P^ 



The next subsection is devoted to the proof of Theorem 14.31 

4.1.1 Proof of Theorem of existence of a solution to the problem (P^) 

Discrete problem. We proceed by using the standard Faedo-Galerkin procedure. 
First of all, let's introduce two increasing sequences Vn and Wn of finite dimensional 
subspaces of V, such that 



|JK = ^ and \JWn = V; (4.10) 



n=0 n=0 



in particular, we choose these subspaces in such a way that Vn C D{B;H) and 

Moreover, we approximate the initial data by appropriate sequences of data such 
that 



n— >+oo 



•(^Cn e Vn Vn, -i^cn > '&0s m ^ 

T1 ^ -j- OO 

XO,n ^Wn Vn, XO,n > Xo lH V 
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The discrete problem (P^) is the following one 

(dtiLogMt) + Xn{t)),v)H + {BMt),v)H = {w{t),v)H 

for a.e. te (0,T), Vt;G K (4.12) 

fi{dtXn{t),u)H + (VXn(t), Vm)^ + (/3^(Xn(t)) + Cx' {Xn{t)) , u) H = {M^) , u) H 

for a.e. t e (0, T), Vm G Vr„ (4.13) 

^n(O) = ?9o,n and Xn(0) = Xo,n a.e. in il. (4.14) 

Theorem 4.4. Leit n ^ N. Under the hypotheses of Theorem \4.3\ and assuming 
^o,n G Vn and xo,n G Wn, the discrete problem (P„) ^J^-^J^ admits a unique 
pair {dni Xn) os solution, such that 

^„eCi([0,r];K) and Xn ^ C\[Q,T]-Wn). (4.15) 

In order to prove this theorem, let {ej}"^! and {6j}"^i be two bases for Vn and 
Wn respectively. Since we can express the functions 'dn and Xn as linear combinations 
of these bases, the true unknowns are the coefficients Uj and yj of such representa- 
tions. If u and y are the vectors which collect these coefficients, then the system 
fl4.12p - fl4.13p may be rewritten in the form of a system of integro-differential ordinary 
equations 

E(t,u(t),y(t),u'(t),y'(t)) = tG(0,T), (4.16) 

where E = (F, G), with components Fj and Gj of F and G defined as follows 

Pj = j Log; ( ^ Uid J ( ^ u'^Ci j(^j + J lYl y'i^i ) ^^ + 5Z ^^ j ^^i ■ ^^^ 

+ a I y^UiCi \ Cj — / wcj (4.17) 

for j = 1, ■ ■ ■ , n; we can notice that the variables (u', y') are independent, as well as 
(u,y). 

The first target is to put the system in a normal form. After that, we will apply 
the Implicit Function Theorem to prove the existence of a solution (at least local). 
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Gj = fj, 



We can notice that E is a continuous function and has continuous derivatives with 
respect to u' and y'. Next, we study its Jacobian matrix with respect to the above 
variables: we can think of it as a four-block matrix, namely the derivatives of F and 
G with respect to u' and y'. Since G does not depend on u', the determinant of the 
Jacobian is 



dE 


= det 


'dF' 
du' 


■det 


'dG' 
_dy'_ 


[d{n',y')\ 



det 



Moreover, all the partial derivatives involved are scalar products: indeed. 



(4.19) 



dy'i 



fx 



bibj 



fi{bi,b 



'j)H, 



dF^ 

du', 




Uk{t)ek CiCj 



^ii ^j)t,\n 



(4.20) 

where (■, ■)t,M is an equivalent scalar product in H , thanks to Lemma WA\ Hence, the 
above matrices are positive definite. 

Now, we find a point (t*, u*,y^, <,y'J such that E(t^, u.,, y^, <, y'J = 0. First 
of all, we set t* = 0; then, we choose u^, and y^ equal to the vectors of coefficients of 
the initial data -(^o.n ^ Vn Xo,n ^ Wn with respect to the chosen bases. 

We define x* n ^s the solution of the following equation 

Axo,n + /3^(xo,n) + cr'(xo,n) = %,n. (4.21) 



/^Xo. 



By comparison, we have x'n ^ -^5 then, we take its projection x**^ on Wn with 
respect to the scalar product in H and we set y'^ equal to the coefficients of Xo*n with 
respect to the basis {bj}. 



Next, in order to find u'^, we define 



'^n{t) = ^Wi(t)ei 



(4.22) 



i=l 



with Wi{t) = {w{t),ei)H] by construction, w„ G C°([0,T];\4) and Wn ^ w in 
C°([0,T]; V"'), a.s n /^ +00. Then, we define «'„ as the solution of the following 
equation 

Log',{'do,n)ul^ + X':n + Bdo.n = Wn{Q). (4.23) 

By comparison, we have u'^^ & H; then, as before, we take its projection on Vn with 
respect to the scalar product (■, ■)t^,u* and we set u'^ equal to the coefficients of the 
projection with respect to the base {cj}. 

Therefore, we have E (t^, u,,, y^, u'^, y'^) = and we can apply the Implicit Func- 
tion Theorem to conclude that the system is (locally) equivalent to a system of the 
form 

(u'(t),y'(t)) = ^(t,u(t),y(t)), (4.24) 
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where <^ is a Lipschitz continuous function with respect to the variables u and y. 
At this point, we can integrate the expression above and we get 

(u(t), y(t)) = (uo,„, yo,J + / ^(s, u(s), y{s))ds =: F{t, u(t), y(t)) (4.25) 

Jo 

where F{t, u(t), y(t)) turns out to be a contraction operator in the space C°([0, r]; R^") 
for < r < T small enough. Using the Contraction Theorem, we claim that there 
exists a unique local solution of the discrete problem. 

Moreover, the solution is indeed a global solution defined on the whole interval 
[0,T]. This result can be proved noting that solution {'&n,Xn) are continuous with 
respect to the time variable; thus, we can consider a new Cauchy problem with initial 
data ('i9„(r), Xn(T)) G Ki x Wn and apply the same method as above to find a new 
pair of solutions with the same regularity of the previous pair which extends it to 
the interval (r, 2r). 

In general, repeating iteratively this procedure, we can find a global solution 

i^n e C\[0,T];Vn) and Xn e C\[0,T];Wn). (4.26) 

Remark 4.5. We could have weakened the hypotheses on the source data and con- 
sidered an inner heat supply with the following regularity g G L^(0,T; H), while the 
regularity of the h-term remains the same. 

In this case, in order to prove the theorem regarding the discrete problem, we 
should have approximated the source operator w G L'^{0,T;V') with a sequence of 
functions {wk} C C°([0,T]; V). 

Furthermore, we identify a new sequence of functions {wn} C C^([0,T];Vn) de- 
fined as 

n 

Wn{t) = '^{wk{t),ei)Hei, (4.27) 

1=1 

so that we can find the point (t*, u^.,y^^, u'^, i/J such that E{t^, u^:,y^,u[,^J = 0. 
For the sake of simplicity, we preferred to introduce a stronger assumption on w. 

Uniform estimates with respect to n. In order to prove the existence of a 
solution to the problem (Pg), e > fixed, we perform some a priori estimates and 
then we let n tend to infinity. 

In this section and in the following ones, we assume the arrangement that the 
symbol c denotes a constant which depends only on the given data; the exact value 
of the constant c may vary in the different estimates and even in the same chain 



of inequalities. A notation like c^ stands for a constant which depends not only on 
the data, but also on the parameter e; nevertheless, it does not depend on the index 
n. We will also use constants like cs, which denotes a dependence on a positive 
parameter 6, but still independent of n. 

The basic estimate is obtained as follows: we test the first equation fl4.12p by 
'(9„(t) and the second equation ( J4.13P by dtXnit), we integrate over (0,t), t G (0,T], 
and we add the equations to each other 



/ dt{Log^d,)^n+ [ |V^„|'+ [ a\^nf + fi [ \di 

jQt JQt J^t JQt 



^tXnf + l I |VXn(t)|' 



n 



+ / PeiXnit)) 



wdr 



(y\Xn)dtXn + / Pe{XO,n) + ^ / |VXo,r 



(4.28) 



For the sake of clearness, we deal with each term separately. First of all, we notice 
that there are some boundary terms (due to the Robin condition that we imposed 
for d) that have to be properly estimated: to do so, we will make use of the Trace 
Theorem. 

Let's start with the first term, which can be rewritten as 

/ dt{h0g,dn)^n= [ leiMt))- //.(^0,n); (4.29) 

JQt Jn Jii 

moreover, the last term is moved on the right hand side of the equality (14.281) and it 
is uniformly bounded, thanks to Lemma [4.21 and (14.111) . 

Thanks to (13. 3p . the second and third terms are estimated from below as follows 



V^„|^+/ al'i?^!^ > c||^„||5^2(oi.y)• 



(4.30) 



The source term on the right hand side of (I4.28P can be easily handled, using 
Holder and Young Inequalities together with the Trace Theorem: 



(4.31) 



g^n+ Mn<5\\^n\\LHQ,fy)+C& 11^11x2(00 + 11^11^^(2.) 



Moreover, we estimate the /3e-term in the following way 



/3.(X0,n) + 77 / |VX0 



Xo,i 



n Jn uq 

< - / |XO,n|^ + - I |VXO,n|^ < C, 

^ Jn ^ Jo 



|VX0,n| 



(4.32) 



19 



thanks to the Lipschitz property of /3^ and (I4.1ip . 

To estimate the a'-term, we have to add on both sides of the equation fl4.28p the 
term 1/2 llXnl^)!!//: on the left hand side a y-norm of Xnit) appears, while on the 
right hand side we perform a chain of inequalities, using (13. 7p . (14. lip and Young 
Inequality, 



1 2 

n\\Xn{t)\\H 



0-'iXn)dtXn 



Xn {dtXn) + 77 IIXO,n||^ " / (r'{Xn)dtXn 

^ jQt 



1 „2 

^ ^ \\XO,n 



H 



Xn {dtXn) + C^ 



Jn 



;i+ \Xn\)dtXn 



Jn 



< C llXo.nlln + ^ \\9tXn\\\2^QA +Cs 

Jo 

< C WXoAh + ^ \\9tXn\\l2,QA +CS WXni 

Jo 



\Xn{s)\\%ds 



H ' 

sjWyds, 



(4.33) 



for each 6 > 0. 

Finally, collecting all the previous results and choosing S > small enough, we 
can apply the Gronwall Lemma (see Reference [161 Appendix, Lemma A. 4]) and get 



Ir/QXIl lIlQll^ ill 11^ ill "^ 

|-'el''^nJllL°o(0,T;Li(f7)) + IP" IIl2(0,T;V') + \\Xn\\ H^(O.T:H) ' \\X' 



^n\\m{0,T;H) 
||Pe(Xn)llL°°(0,T;Li(n)) — '^e- 



n\\L°°(0,T;V) 



(4.34) 



after having taken the supremum over t G (0,T). 

Next, we derive an a priori bound for dt^n ^ -^^(0, T; H); we test (I4.12p by dt^nit) 
and integrate over (0,t): 

L0g',{^n)mnf+ I {dtXn)mn) + \\\Vd^{t)\\\+ [ a^df^n) 
jQt ^ J^t 

= f 9ndt^n+ f hndtdn + l\\V\n\\l. (4.35) 

jQt JT,t ^ 



Thanks to Lemma 14.11 the first term is easily estimated 

Log^(^n) \dt^n\^ > e ||<9i^n|lL2(0,t;i?) • 

Concerning the second term, we move it the right hand side and we get 



(4.36) 



dtXn dt^n 
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(4.37) 



where we applied Holder and Young Inequalities; moreover, the norm of dtXn is 
uniformly bounded, thanks to the previous estimate (14.341) . 

Due to (13. 3p . (14. lip and the Trace Theorem, the boundary term is treated as 
follows 

aMdt^n) > I WUml^r) - c \\A,n\\l > ? WUml^iD - c. (4.38) 

Finally, after having noticed that || Vt^o.nllj:^ is uniformly bounded (see (14. lip ), we 
estimate the source terms: 



/ gdt'&n < ^\\9tl3n\\L'2{0,t;H)+ '^s\\9\\l2{0^T;H)'-> 
J Of 



hdt^n = / h{t)Mt) - / HO)^o,n - / {dth)^n 

St 7r Jr J'Et 

\dMs)\\^,^^^)\\Us)\\mT)ds; (4.39) 

we recall that h G L'^i'd.T- L'^iV)) nW^^^iQ.T; L'^iV)). 

Collecting all the estimates above, we can choose 5 > Q small enough and apply 
the Gronwall Lemma (see Reference ^^ Appendix, Lemma A.5]), in order to get 

ll^n|lL-(0,r;y) + \Wn\\l2(Q^T;H) < ^s- (4.40) 

Passage to the limit as n /^ +00. Using the previous estimates, we can claim 
that there exist 

^e^H\id,T-H)fM'^{Q,T-V) and Xe e i7^(0,T; i7) n L°°(0, T; \/) (4.41) 

such that ^n -^ ^e and Xn -^ Xe weakly in if^(0,T; H) n L°°(0,T; V), as n / +cx) 
(at least for a subsequence). 

These weak convergences are already sufficient to ensure that the Cauchy condi- 
tions hold in the limit: '(9e(0) = T^o.e and Xe(0) = Xo,At i^ H and a.e. in ^2. Moreover, 
applying Aubin Lemma (see Reference [161 Appendix, Lemma A.5]), we have 

^n^^e and Xn^Xe mL\0,T;H). (4.42) 

Strong convergence implies pointwise convergence (a.e.) in Q; therefore, we can 
identify the limits of the non-linear terms, taking into account the Lipschitz property 
of /3'^, Log^ and a'. 
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At this point, we can show that the pair {i)^, Xe) is indeed a solution to problem 
(Pe), e > 0. Fix n G N and test equations fl4.12p and fl4.13p by two arbitrary functions 
V G L^(0,T; Vm) and u G L^(0,T; Wm) respectively, with m < n; then, we integrate 
over (0,T). Now, we let n tend to infinity and, using the proved convergences, we 
get 



dt(Log^'&e + Xe)v + / V-i^e ■ Vf + / a-i^^f = gv + hv 

-^^ ^^ ^"^ ^) (4.43) 

^Qt Jq Jq Jq Jq 

Wv G L^(0,T; \/m) and Vit G L^(0,T; 1^™.), Vm G N. Since m is arbitrary, the same 
variational equations hold for each v G L^(0,T; V), by density. 

Moreover, by comparison in the second equation of ( I4.43p . we have —Axe ^ 
L'^{0,T; H), then, due to Elliptic Regularity Theorem (see [HI Chap. 1.5]), we can 
conclude 

XeeL\0,T;DiA;H)). (4.44) 

Uniqueness of solution to the problem (P^). In order to prove the continuous 
dependence of solutions from data (hence, uniqueness), we follow the same method 
illustrated in [I] Sec. 5]. 

Let consider the time- integrated version of (14.60 . namely 

Log^C + Xe + l*5^£ = l*w + '7o,e where ?7o,e := Log^-i^ce + Xo,M5 (4-45) 

and we couple it with the second equation (14.71) of (P^). We pick two solutions 
{^s,i,Xe,i) to the system corresponding to the sets of data (w7j,?7o,£,i), i = 1,2. Then 
we write both (I4.45P and (14. 7p for such solutions and multiply the difference of 
the first equations by -i^e := i)^,! ~ ^e,2 and the difference of the second ones by 
Xe '■= Xe,i ~ Xe,2- Finally, we sum the obtained equalities to each other and integrate 
over Qt := Q X (0,t). 

After some manipulations, we have 

/ [LogAi-LogA,2R + ^l|l*V^,(t)||?,+ / (l*«C)C + ^||x.(t)||?, 



+ / |VXe|'+ / [/3;(x.,i)-/3^(Xe,2)]x.= / (l*^ + r7o,.)^. 

jQt JQt JQt 

- [ [cr'{Xe,i) - cx'{xe,2)]Xe + ^ Wxojl ' (4-46) 

JQt ^ 



'Qt 

'Qt 
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where we introduced a similar notation for all the differences involved: w := wi —W2, 

XO,fi '■= XO,fJ.l — X0,At2- 

We deal with each term separately. Concerning the non linear terms, we have that 
the Log^-term and the /9^-term are non-negative, since Log^ and /3^ are monotone, 
while the a'-term is estimated as follow (see (13.611 ) 



(4.47) 



/ W'iXe,l) -(^'iXe,2)]Xe<CL \Xef = Cl \\Xeis)\\% ds. 

JQt JQt JO 

The boundary term can be easily handled 

[l*a^e)^,>^\\l*Mt)\\Ur)- 
Regarding the first term on the right hand side of the equality, we have 



(4.48) 



il*W + r]o,e)^e= / il*9it)+VoJil*Mt))+ / (l*/i(t))(l*^e(t)) 

Jn Jv 

- [ g{l*^e)- [ h{l*^,)<cs(\\l*gml + \\mAl+U*Ht)\\l,^^^ 

JQt JT.t ^ ' 

+(5||l*^,(t)||?,+ /" (||^(5)||H+IIM^)llL^(r))l|l*^.(^)lly, (4.49) 



for all (5 > 0, thanks to Young and Holder inequalities and Trace Theorem; we recall 
that g e L'^{0,T; H) and h G L^(0,T; L^(r)). Moreover, using Young Theorem (see 
Reference [261 Chap. 8, Prop. 8.9]) with r = oo and p = q = 2, we have 



111 * 9{t)\\H + \\Vo,efH + P * Ht)\\l2(^r) < 111 * ^llL-(o,r;H) + WvoAh 

|1 * ^llL°°(0.T;L2fr)l - '^ 



I l|2 I II ||2 I ||t,||2 

\9\\l2(0,T;H) + \\VO,e\\H + ll"'llL2(o,T;L2(r)) 



(4.50) 



Now, we choose 6 > small enough, so that we can apply Gronwall Lemma (see 
Reference [Ml Appendix, Lemma A. 5]) 



f [LogA,l-Log,^,,2]^.+ /" mXe,l)-PeiXe,2)]Xe+\\VXe\\l2 
JQt JQt 



(Qt) 



+f^\\Xeit)\\%+\\l*l3, 



V 



< c 



Q iiii9 iiii9 

l'70,e||j/ + \\9\\l2{0,T;H) + II^IIl2(0,T: 



L2(r)) 



(4.51) 



Vt G (0, T). This concludes the proof of Theorem 14. 3 [ 
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4.1.2 Convergence of problems (P^) to the original problem (P^) as £ \ 

Theorem 14.31 ensures that for e > fixed there exists a solution to the approximated 
problem (P^) and it is unique. Now, we want to get a solution to the original problem 
(P^) letting e tend to zero. We consider an arbitrary solution ('i?^, Xe) of (Pe) and we 
will perform a number of a priori estimates so that we can take the limit as e \ 0. In 
general, these estimates will hold for e small enough; however, they are independent 
on the parameter e. 

First a priori estimate. The first estimate is analogous to the first estimate 
performed in the calculations above. As before, we test the first equation (14.61) by 
^eit) and the second one (14. 7p by dtXeif)] then, we integrate over the time interval 
(0, t) and we sum the obtained equalities to each other. 

We can notice that all the terms are bounded by a constant c, independent of e, 
except for the /S^-term, which is estimated as follow 

/3.(Xo,m) < / /3(Xo,m) < c, (4.52) 

due to the properties of Yosida regularization and to the assumption (13.81) . 
Then, we have 

+ ll/^£(X.)llLoo(o,T;Li(n)) < c. (4.53) 

and, by comparison, \\dtLog^{'de)\\i2iQT-v') — ^■ 

Second a priori estimate. The following estimate provides us a bound for the 
non linear term (3'^. We test the equation ( 14. 7p by [i'eiXe) and we integrate over (0, t), 

te(o,T] 

l,JMXe{t))+ [ VXs-VmXs)]+ [ mXe)f= [ ^eP'eiXe) 

Jn JQt JQt JQt 

- I (y'{Xe)Pe{Xe) + /i / /^.(Xo,^)- (4-54) 

JQt Jn 

The second term on the left hand side is non-negative, since /3^ is monotone; the 
right hand side is estimated in the following way 

^el3'e{Xe)- [ a' (Xe) (3'eiXe) + l^ [ MXo) 

JQt Jn 

< ^ WeiXs)\\l2^o^f^H) + C5 (||^.|li2(o,r;H) + \\Xe\\l2^o^T;H)) + ^ (4.55) 
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for all 6 > 0, due to Holder and Young Inequalities, (13.6^ . (I3.8P and < /3e(r) < /3(r), 
V r G R. We recall that the norms of -^s and Xs are uniformly bounded, thanks to 

Then, we choose 5 > small enough and we take the supremum over t G (0,T): 



\\l3eiXe)\\l2(^0,T;H) < C 

Collecting all the estimates founded so far, we have the estimate 

\\Xs\\l2{0,T;D(A;H)) — ^■ 



(4.56) 



(4.57) 



Third a priori estimate. The last estimate we need is a bound for the logarithmic 
term in a suitable functional space. We recall that, by comparison in the first equa- 
tion (14. 6p . we already know that Log^('i9e) G if^(0,T; \^'), with uniformly bounded 
norm (with respect to e) thanks to (I4.53p . 

In order to obtain this estimate, we have to be very careful in dealing with the 
boundary terms, so that we can obtain a suitable bound for the logarithm. 

Remark 4.6. We recall that, since Log^{r) /^ log(r), Vr > 0, as e — t- 0+ (see 
Reference f7^ Chap. II]), then 



Log^{r) < r, Vr > 1 and Log^{r) < 0, Vr < 1. 



(4.5^ 



We test (14.61) by hog^i^e) and we integrate over (0,t) 

^||Log,^,(t)l|^+ / V^,-V(Log,^,)+ f «^,(Log,^,) 
= [ g{Log,i^,)+ [ /i(LogA)- / {dtXe)l^og,^, + l\\Log,do,e\\l. (4.59) 

jQt JT,t JQt ^ 

As Logg is a monotone function, the second term is non negative; the boundary 
term needs a deeper analysis: keeping into account Remark 14.61 we have 



St ^0 



L^rn{i?ee(-oo,o)u(i,+oo)} "'rn{i?ee[o,i]} 

t 



>a 



^sLog^^, 



rn{i9,6[o,i]} 



c'^e 



> -c, 



(4.60) 
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where the last inequahty follows from a suitable application of Lebesgue Convergence 
Theorem. 

Concerning the right hand side of the equality, the source terms are treated in 
the following way: the g-teim can be easily handled together with the 9tXe-term 

/ ^(LogA)- / (9,x.)LogA< / i\\9is)\\H+\\dtXeis)\\H)\\l^ogMs)\\H^ 
JQt JQt Jo 

(4.61) 

we recall that g G i^^(0, T; H) and dtXe ^ -^^(0, T; iJ), with uniformly bounded norm, 

thanks to f l4.53p . On the other hand, the /i-term requires some additional calculus 

/ /i(Log,0</ / ^^£ < ll'^e|li2(o,T-L2(r)) + ll^llL2(o,T-L2(r)) <c, (4.62) 

recalling that h> a.e. on S and using Trace Theorem and (I4.53p . 

Regarding the last term, provided that e is small enough (say e G (0, 1)), we have 
the following bound 

l|Loge'^o,£|lioo(j^) < c II VllL-(n) ^ c. (4.63) 

using properties stated in Remark 14.61 and thanks to (14. 5p . 

At this point, we can apply Gronwall Lemma (see Reference [T6l Appendix, 
Lemma A. 5]) and take the supremum over t G (0,T) 

\\Log^l^e\\L--(0,T;H) < C. (4.64) 

Passage to the limit as e \ 0. Collecting all the previous estimates and using 
suitable compactness results, we can state that there exist 

X, G L'iO, T; D{A- H)) n H\0, T; H) 

^ G L°°(0, T; H) D H\0, T; V) 

such that, at least for a subsequence e„ \ 0, they are limit of the approximating 
solutions 

^, -^^ ■mL^iO,T;V) (4.66) 

Xe-XM inL\0,T;D{A;H))nH\0,T;H) (4.67) 

/3:{Xe)^^^. mL\0,T;H) (4.68) 

Log,(C)^^ mL°°{0,T;H)nH\0,T;V'). (4.69) 
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It is clear that the Cauchy conditions (J3.12P are fulfilled. Indeed, weak conver- 
gence in H^{0,T]X) of a sequence implies weak convergence in X (any functional 
space) of the corresponding initial data. Thus we can directly conclude Xm(0) = Xo,^ 
in H and a.e. in Q. 

For the logarithmic term, we have Logg'i9e(0) = Log^-i^o.e ^ -^(0) i^i ^'; on the 
other hand, Log^'i^o.e -^ logf^o for a.e. x & Q (hence, it converges in measure) and 
Log^t^o.e £ -^°°(^); with uniformly bounded norm. Then, Vg G [1, +00) we have (see 
Reference [Sni Chap. XI, Proposition 3.10]) 

LogA,. ^log^o inL\n); (4.70) 

in particular, Log^-i^o.e "^ log-f^o in L'^{Q). Then, by uniqueness of limit, -2'(0) = 
log-i^o in H and a.e. in Q. 

Now, we deal with the non linear terms. Thanks to Aubin Lemma (see |T71 Chap. 
1.5, Theorem 5.1]), Xe -^ Xm in L^{0,T;V), as e \ 0, then a'{xe) -^ cr'iXp.) in 
L'^{0,T] H), since a' is Lipschitz. As far as the non linearity associated with the 
maximal monotone graph (9/3 is concerned, we can apply Barbu Lemma (see [T5| 
Chap. II, Lemma 1.3]) and conclude that x^l ^ D{df3) and ^^ G dPixti) a.e. in Q. 
Finally, concerning the logarithmic term, thanks to the compact embedding if C V^' 
and Aubin Lemma (see [171 Chap. 1.5, Theorem 5.1]), we have 



l^og^^e -^ ^ *-weakly in L°°(0, T; E) and strongly in ^^(0, T; V% (4.71) 

thus, log^(?9£) -^ S£ weakly in L^{0,T]H) and strongly in L\0,T;V'), as £ \ 0, 
since ei)s ~^ strongly in L^(0,T; V). Therefore, the following limit holds 

hm/ {\ogM),^e)=f i^,^,); (4.72) 

^^"Jo Jo 

where (■, ■) stands for the duality pairing between V and V; we apply Barbu Lemma 
once more and we conclude {}^ G -D(log) (then, {}^ > 0) and ^ = log -(9^ a.e. in Q. 

Remark 4.7. We can notice that we could have taken two different parameters e 
and e' in approximating log and dfi, in order to prove Theorem \3.4\ 

Moreover, we could have kept fixed either of them, say e, and let e' tend to zero. 
This leads to an existence result for a semi-regularized problem. All the a priori 
estimates are conserved in the limit (thanks to the semicontinuity of the norms), 
so that it is possible to let e tend to zero. The same can be done by switching the 
parameters. 
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Uniqueness of solution to the problem (P^) and continuous dependence on 
data. We conclude the proof of Theorem l3.4l by showing uniqueness and continuous 
dependence of solutions from data. We follow the same procedure used previously 
to prove uniqueness of solution to problem (Pe). 

Let consider the first equation fl3.9p . formally integrated in time, and the second 
equation ( ]3.10p of problem (P^) for two solutions {'&^,i, Xii,i: ^^l,i) corresponding to the 
sets of data {wi, ^o,i, Xo,/.,j), i = 1,2. 

We test the difference of the first equations by -i?^ := "i?^,! —'&^,2 and the difference 
of the second ones by x^l '■= Xfi,i ~ Xn,2] then, we integrate over (0,t) and we sum 
the obtained equalities 

/ [log^^,i-log^^,2]^^+ / (1*V^^)V^^+ / (l*a^^)^^ + /i [ {dtx,)x, 

jQt JQt •iT.t JQt 

+ / \^xX+ I i^.X,= I (l*w^ + ^oJ^M- / [^'(Xm)-^'(Xm,2)]X/.. (4.73) 
JQt JQt JQt JQt 

where we set again a similar notation for all the differences involved: w := wi — W2, 
After appropriate estimates, we get the following 

/ [logV-log^M,2]^M+ / ^^^X^ + \\'^xX^iQt)+^'\\x^^(^)\\l + \\'^*Mt)\\l 
JQt JQt 

< M [\\vqJ\h + \\9\\l^(o,T;H) + ll^llL2(o,T;L2(r))) ' (4-74) 

with M = M{Q, T) and g := gi — g2, h := hi — /i2- 

From this relation, continuous dependence on data follows immediately and unique- 
ness of solution too (set gi = g2, hi = /i2 and r/c/x,! = Vo,fj.,2)- 

5 Asymptotic behavior as /i \ 

In this section we prove Theorem 13.61 and study the asymptotic behavior of the 
problem (P^) as the time relaxation parameter /x \ 0. 

The procedure is the following. We first prove some a priori estimates, uniform 
with respect to /i, and then let /x tends to zero. Using suitable compactness results, 
we find a subsequence of solutions to problem (P^) which (weakly) converges to a 
solution to the limit problem (Pq); the whole family of solutions indeed converges, 
thanks to the uniqueness result stated in Theorem 13. 7[ which is proved at the end of 
this section. 
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We make a distinction between the two cases (Hypothesis 1 or 2), since the 
passage to the hmit as /i \ and the proof of existence of solutions to the hmit 
problem use different methods in either case. 

First of all, let point out that all the a priori estimates, which we perform, are 
formal: we should come back to the procedure used in pLj and in the previous section, 
where problem (P^), /i > 0, has been solved by letting e \ in the approximating 
problem (Pe). However, in order not to make the exposition too heavy, we prefer to 
proceed formally. 

All the constants that appear in the following estimates are positive quantities, 
independent of the parameter /i. 

5.1 Uniform estimates 
5.1.1 First a priori estimate 

The first estimate gives a uniform bound for the corresponding norms of solutions 

Assume either Hypothesis 1 or 2. We formally multiply equation (13. 9 p by ^^ and 
equation (13.101) by dtXfj., then we integrate over (0,t), with t G (0,T], and we sum 
the resulting equations. After some rearrangements, we have 

Jn JQt JT.t JQt ^ Jn JQt 

+ [ cr\x,)dtX, = h,+ I w^, + ]- I |Vxo,mI'- (5.1) 

JQt Jn JQt ^ Jn 



Let deal with each term separately. The first term on the left hand side is non 
negative, since d^ > Q a.e. on Q; the second and third terms are estimated from 
below by the norm of d^ in L^(0, t; V). 

For the last two terms of the left hand side, we distinguish the cases in which 
Hypothesis 1 or 2 holds. Under Hypothesis 1, using a' = a constant and (I3.27p . we 
have 

^Ax^l + / (^' {Xt^)dtXfi = / f^iXf^ii)) + « / X^W " / /3(Xo,m) - a / Xo,/. 
t JQt Jn Jn Jn Jn 

> ci \\Xfi{t)\\\ - C2 + a / x^l{t) - / /3(Xo,m) - « / Xo,m; (5-2) 

Jn Jn Jn 
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we move the last four terms on the right hand side of the inequahty (I5.ip and, thanks 
to fl3.26p and Young Inequahty, we estimate them this way 



(5.3) 



-a X,.{t) + / /3(xo,^) + a I xo,m < S Hx^WII?/ + cs V5 > 0. 
Under Hypothesis 2, using fl3.29p . we have 

^AX^. + I (y'{x^.)^tX^. = I l^ixM + / <xM - [ /3(xo,m) 
jQt Jn Jn Jn 

- / (^(Xo,m) > ci 11x^(^)11?/ - C2 - / /3(xo,m) - / ^(Xo,m), (5.4) 

Jn Jn Jn 

where the last two terms are uniformly bounded, thanks to ( 13. 7p and ( 13.26^ . 

Regarding the right hand side of equality (15.11) . the last term is uniformly bounded, 
thanks to ( 13.26^ . while the w-term is bounded from above by the following norms 

(5.5) 

hl^^, < 6 11^^11^2(2,) + C5 ||/i|lL2(2,) < S ||^/.||i2(o,i;y) + Q 11/^11^2(2,) , 

using Holder and Young Inequalities and the Trace Theorem. 

Then, choosing 6 > sufficiently small and taking the supremum on [0,T], we 
conclude 

\\^tJ-\\L'^{0,T;V) "*" ^^ ll^*^MllL2(o,r;_H') "*" H^M llL°=(0,T;y) — ^' \^-^) 

moreover, by comparison, we have ||log^„|Loo/nTvo — ^■ 



A'HL°°(0,r;V) 



5.1.2 Second a priori estimate 



In order to obtain a uniform estimate for the norm of the non linear term ^^ G dP{x^), 
let formally multiply the equation (I3.10p by ^^ and integrate it over (0,t), t G (0,T]; 
after some rearrangements, we get 



/i / /3(XmW)+ / (3"iX 



n JQt JQt Jn jQt JQt 

(5.7) 
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The first term on the right hand side of this equahty is uniformly bounded, thanks 
to (13.261) . while we can easily handle the other two terms applying Holder and Young 
Inequahties 



Qt 



(5. 



for each arbitrary 5 > 0; we recall that a' is a Lipschitz function (see (13.61) ) and that 
the norms of -(9^ and Xn ^-^e uniformly bounded with respect to n thanks to estimate 

We choose 6 > small enough and we take the supremum over t G (0, T); finally, 
we have in particular 

moreover, by comparison in the equation (I3.10p and by Elliptic Regularity Theorem, 
we can conclude ||XmIIl2(o,t;D{A;//)) ^ c. 

5.1.3 Cauchy estimate 

Consider two problems (P^) and {Pi,), with /i > z/ > 0; repeating the same arguments 
as above to show the uniqueness of solutions to problem (P^) ( Section H]), subtract 
the first two equations (13. 9 p integrated in time and the second two equations (I3.10p . 
Let multiply the resulting expression by -(^^j — i!}^, and by x^ — Xu respectively; 
finally, we integrate over the interval (0, t) and we sum the equations 

/ [log^^-log^,](^^-^,)+ / (1*V(^^-^.))V(^^-^.) 

JQt JQt 



+ {l*a{^,~A)){^,-A) + X dt{x^-Xu){x,.-X.)+ |V(XM-X.)r 
JT,t JQt JQt 



+ 1 {^i^-^y){Xt,-Xv)+ j [(y'{Xt,)-(^'{Xy)\{Xii-Xv) 

JQt 

< I (Xo,M-Xo,.)(^M-^-)' (5-10) 



for some z/ < A < /i. 

At this point, we distinguish the case where Hypothesis 1 or 2 holds. 
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Under Hj^othesis 1, a' = a G R constant, which imphes that the term involving 
o' is identically zero. 

After some manipulations and estimates, we get the following inequality 



[log ^, - log ^ J (^^ - ^,) + / (e^ - i^) {x, - Xu) 

jQt 

+ 111 * (^M - ^.')Wlly + II VXm - Vx.||i2(Q,) < c ||xo,M - Xo,.|Ih , (5-11) 

vte (o,r). 

Thanks to this estimate and to (13.261) . we can conclude that {1 * t9^} and {Vx^} 
are Cauchy sequence in L°°(0,T; V) and L^(0,T; H) respectively. 

Under Hypothesis 2, using the strong monotonicity property f l3.28p assumed for 
(9/3 + 0"', we have 



JQt JQt JQt 

> I|V(Xm - X^)lli2(o,T;H) + P WXpt - X^ 1 1^2(0,7;//) > c llx/. - xA?L2^Q^T■y) ■ (5-12) 

Then, after some estimates, we get 

/ [log^^ - log^^](^^ - ^^) + 111 * (^^ - "&u)\\\^(Q^T;V) + IIXm - Xu\\\2(^o,T;V) 

J Qt 

< c||xo,m-Xo,!.|Ih- (5-13) 



and we conclude that {1 * 'd^^ and {x^} are Cauchy sequence in L°°(0,T;y) and 
L^(0, T; V^) respectively. 

5.2 Passage to the limit as /i \ under Hypothesis 1 

Using the previous estimates, we can state that there exist 

i^L\<d,T-H) ^ ' 

CgL~(0,T;F') 
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such that the following convergences hold (at least for a subsequence) 

■(9^ ^ ^ in L\0,T;V) (5.15) 

1*^^ -^ 1*^ inL°°(0,T;1/) (5.16) 

X, ^ X inL\0,T;D{A;H)) (5.17) 

MXm ^ mL\0,T;H) (5.18) 

^, - e mL\0,T;H) (5.19) 

log(t9^) - C mL^iO,T;V'). (5.20) 

Now, we deal with the non linear terms. We mainly refer to [21 Sec. 3.3], where 
an asymptotic behavior as the interfacial energy coefficient tends to is studied for 
the same PDE system with different boundary conditions. 

In particular, in order to analyze the term with the maximal monotone graph (9/3, 
we perform an ad hoc weighted estimate, choosing t'^'^ as weight, with 7 = 3/4. The 
choose of such a weight will be clear in some later inequalities (see (I5.34p ). which 
happen to be fulffiled if 1/2 < 7 < 1, and we take 7 = 3/4 at once. 

Lemma 5.1 (see Lemma 3.1, [2]). Let 7 > 0, then the following uniform estimate 
holds 



^2,^" 



hL+ sup t'-^w^.ml + fi sup t^^ f \d,xM'+ [ t'^ \vda,f < c. 

*'/u te(o,T) te(o,r) Jn Jq 

(5.21) 



Proof. We formally test the equation (13. 9p by t^'^dt'dfj,] next, we formally differentiate 
(I3.10p with respect to time (we obtain a second order equation) and we test it by 

We notice that, since a' = a, then a" = 0. Finally, we add the equalities to each 
other and we integrate over (0,t): 

dtlog^,t^^dt^,+ [ V^,t'^V{dt^,)+ f a^,t^^dt^, + fi [ t^WtXf^dh, 
jQt J^t JQt 

+ I t'^\y{dtx,)?+ I t'^i',\dtxX= I gt'^dtd,+ [ ht'Wt^,. (5.22) 
JQt JQt JQt J^t 

We deal with each term separately. 

The ffist and the last two terms are non negative, since (recalling that i}^ > 0) 

/ {dtlog^,)t'^{dt^,)= [ t27M^>o, (5.23) 

JQt JQt *^M 
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and dfi is nionotone;next, we estimate the second one as follows: thanks to (]5.6p . 

^27 



jQt ^ J^ JQt 

+27 +27 

^ l|V^,(t)||' - 27T2^-^ \\V^A'^.„, > — \mM)\\l - c. 



Similarly, using (15. 6p . we have 



(5.24) 



^27 
>a — M 



^27 



and 



2 „"Ml'';ili2(r) - 2a^T^^ ^ \WJ%{o,T;LHt)) > "^ ll^^^^)ll^^(r) ~ ^' (^-^S) 



+27 



JQt ^ Jn JQt 



/it^^ 



/it^T 



The source terms can be easily handled 

/ gfWt^, = e^ f g{t)d^{t) - f 2^f''-Xg - f e^d.dtg 



<5e'U„{ml + C5\\g{t)\\\ + cT^^-' 



Mi^-JiiH 



9 ||2 I II ||2 

^MllL2(0,T;fl') + \\9\\m{0,T-H) 



<^t^'¥M\l + c, (5.27) 

for each 5 > 0; in the last inequality we used g G H^{0,T; H) C C°([0,T];if) 
and (15.61) . The same calculations hold for the h term, thanks to the fact that h G 
H\0,T;L'^{T)) C C°([0,T]; ^^(r)), the Trace Theorem and (l53|) 



/it'^5i^;.<5t27||^^(t)||;,2,r,+C, 



M'^''^llL2(r) 



(5.28) 



for each 6 > 0. 

Choosing 6 > small enough, we have 

f t'^\^ + t'^f\\VMml + \\MmUr)+f^\\dxMfH 
JQt ^1^ ^ 

+ / t'^\V{dtX,)?<c, (5.29) 

JQt 

and we get the desired result by taking the supremum over t G (0,T). D 
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At this point, we use the identity 
and observe that, thanks to Sobolev Inequahties and (15. 6p . 



(5.30) 



t^M 



L°°(0,T;L12(J7)) 



llr^, 



MIIL°°(0,T;L6(Q)) — '^-^^ \\^fJ-\\L^{0,T;V) — ^■ 



Thanks to the previous Lemma, we can say [[f^ {(^t\f&^ \\L'^inT-m 
look for exponents p,q,r,s > 1 such that 



(5.31) 
< c. Now, we 



t^(5tv^)-t^v^-t" 



< 



^^ {dtV^) 



Li{0,T:L^(n)) 



L'2{0,T;H) 



t^M 



L°°(0,T;L12(Q)) 



t 2 



L=(0,T) 



<C, 



(5.32) 



and 



Iy-^"^i9 II < lly-^^MI Ili9 II 

I'' '^A«||LP(0,T;L6(n)) - ir IIl-CO.T) 1 1*^^11 L2(0,T;Le(Q)) 

< ll/T^^II IU9 II < r 

- Ir IIl'-(o,t) II'^a'IIl2{o,T;V) - ^' 



(5.33) 



where we used Holder and Sobolev Inequalities. Thus, we need the following con- 
straints to be fulfilled: 

- = l + -, - = l + -. t-iGL^(0,T), t^-'eL^{0,T). (5.34) 

q 2 s p 2 r 

As 7 = 3/4, we can take p = 6/5, r = 3, q = 14/13, s = 7/3. Then we have 

\\'^t{t ''^/i)llLi4/i3(o,T;Ll2/7(Q)) < C. (5.35) 

At this point, we can apply a standard compactness lemma (see Reference [H Sec. 
8, Corollary 4]) and obtain tt (^^^ _ ^) ^ o strongly in C^{[0,T]; H); then. 



^mXm^^X weakly in L^(Q 



t), 



(5.36) 



since, ^^x^. = t"^^^t~^'\^., X/. ^ X m i:°°(0, T; r) and r^'^ E L\0, T). 
Finally, collecting all the previous results, we have 



limsup / ^^Xm < / ^X- 
M^o Jq Jq 



(5.37) 
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Indeed, by comparison in the second equation ( I3.10p of problem (P^), we have 
Jo Jo Jo Jo Jo 



'Q JQ JQ 

_ II ||2 _ 11^ 

2 \\XQ,^l\\H II ^X/iiii2(Q) 

||2 
^A'IIl2(Q) 

On the other hand, thanks to (]5.17p . (J5.36P and the hnearity of a 



^ ^ II l|2 IIV7 l|2 / I / Q 

< o IIXcmIIj^ - l|Vx^lli2(Q) - I ax^+ I ^f^X,^ 

</^c- II Vx^ 11^2(0) - / «Xm+ / ^mXm- (5-38) 



Ux = - I I Vxl' - I ax+ hx, (5.39) 

7q jQ Jq Jq 

hence 

hmsup /" e^XM <- I |Vx|' - / ax + / ^X = I ix- (5.40) 

M-^o jQ jQ Jq Jq Jq 

In conclusion, owing to [15j Sec. II, Lemma 1.3], we get x ^ D{df3) and ^ G df3{x) 
a.e. in Q. 

Concerning the logarithmic term, we can follow the same procedure we have seen 
for the ^ term. By comparison in the first equation integrated in time and due to 
the previous convergences (see f IS.lSp . fl5.16p . f l5.36p ). we get 



hmsup [ {C„^^)< [ (C,^), 
/i-i-O Jo Jo 



(5.41) 



where C/^ = log"*?^. Hence, owing again to [151 Sec. II, Lemma 1.3], and knowing that 
C;. e Logt?^ = <9*(^^) (see Remark 4.3]), we have ^ G D(Log) and C G 9^(?9) = 
Log ■(9. 

It is also possible to prove that the solution ^ to problem (Pq), which is the 
absolute temperature, is strictly positive (see pi Theorem 4.7]). 
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5.3 Passage to the limit as ^u \ under Hypothesis 2 

Taking into account the uniform estimates we have performed, we can conclude that 
there exist 

i^eL\0,T;V) 

XeL\0,T;D{A;H)) 

^eL\0,T;H) ^ ■ ^ 

CgL~(0,T;\/') 

such that the following convergences (as /i \ 0) hold, at least for a subsequence, 

1^^ ^ ^ in L'^{0,T;V) (5.43) 

1*^^ -^ 1*^ mL°°{0,T;V) (5.44) 

X^ ^ X weakly in L2(0,T;D(A;i7)) 

strongly in L2(0,T;\/) (5.45) 

fidtx>. ^ mL\0,T;H) (5.46) 

e^ - e mL'iO,T;H) (5.47) 

log(^^) - C inL-(0,T;n- (5-48) 

Thanks to these convergences, we can identify the a' term and we can immediately 
apply the result stated in [151 Sec. II, Lemma 1.3], in order to get x ^ D^dfi) and 
^ G /^(Log), i E d(3{x) and C E 9^(^) = Log^ a.e. in Q. 

5.4 Uniqueness of solution to problem (Pq) and continuous 
dependence on data 

To show the uniqueness of solution and continuous dependence on data, we follow 
the same method used in Section HJ to show uniqueness of solution to the problem 
(P^) with n>0. 

In particular, we distinguish the case where Hypothesis 1 or 2 holds, since the 
behavior of the non linear terms /3 and a is quite different. 

Assuming Hypothesis 1, we write problem (Pq) for two distinct solution {{}i, Xi, C,i, 0); 
i = 1,2, and we subtract one equation to the other. 

[Ci - C2] + [Xi - X2] + 1 * 5(^1 - ^2) = 1 * {wi - W2) + r7o,i - r7o,2 (5.49) 

A{xi - X2) + [^1 - 6] + [^'(Xi) - cr'{x2)] = ^1 - ^2. (5.50) 
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We test ( I5.49P hj {} := {}i — {}2 and (15.50^ by x •= Xi ~ X2 and we integrate over 
(0,t), t G (0,T]; finally we add the resulting equation: 



Vo)^, 



f C^+ [ (1*V^)V?9+ / {l*a^)^+ f \Vxf+ f ^X= f {l*w + 
jQt JQt J^t JQt JQt JQt 

(5.51) 

where ( := (^ - (3, ^ := ^1 - 6, w := Wi - W2, ^0 := ^0,1 - ^0,2 ; the a' term is 
identically zero, since a' = a constant. 

After suitable estimates, similar to those performed in Section IH we get the 
following inequality Vt G (0,T) 



C^+ / ex+l|Vxll2.(Q,) + l|i*^(t)ry 

< ^ (holli + l|t/||i2(o,r;H) + ll^llL2(o,T;L2(r))) > (5-52) 
where M = M(f2, T) and g := gi — g2, h := hi — /12. 

Assuming Hypothesis 2 and using (13.281) . the terms related to d(3 e a are handled 
as follows 

f Nxf+ f ^X+ f W iXi) - ^' iX2)]x >c Ml. ^o,T-y)- (5-53) 

jQt JQt JQt 

For the other terms we apply the same procedure seen in Section H] and we get 

/ C"*^ + ||1 * ^llL°°(0,T;y) + 11X11^2(0,^;^) 

•^ Qt 

< M (\\r]o\\\ + il^lli2(o,T;//) + ll^lli2(o,T;L2(r))) ' (5-54) 

where M = M((],T). 

Finally, we can notice that, since \E' is strictly convex, its subdifferential d'^ is 
strictly monotone, then we have t?i = '§2 and Xi = X2, if we set gi = g2, hi = /i2 and 

^0,1 = '70,2- 
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